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Abstract 

We calculate the one-point function of certain chiral operators, known as open Wilson lines, 
in the holographic Nambu-Jona-Lasinio (NJL) model at finite temperature as well as finite 
chemical potential for quark (baryon) number density. These operators are holographically 
dual to Euclidean string worldsheets bounded by the flavor branes. Their one-point func- 
tions can serve as order parameters for chiral symmetry breaking in these models. Analyzing 
the behavior of the dual string worldsheets with respect to temperature and chemical poten- 
tial (or quark density) enables us to determine how the order parameter depends on these 
parameters. 



1 Introduction and Summary 



QCD at finite quark density and temperature has a far more interesting phase structure than 
at just finite temperature; see [1] for a review. For instance, depending on the number of 
flavors, new phases such as color-superconductivity and color-flavor locked are believed to 
emerge at moderate to high densities. While the phase diagram is well understood at zero 
density and finite temperature, little is known at finite density and temperature, especially 
in the regime of intermediate densities where QCD is strongly coupled. Needless to say, 
understanding the phase structure of QCD at moderate densities will help us understand 
the physics of the quark-gluon plasma produced at RHIC, as well as the physical processes 
within the core of neutron stars. Unfortunately, lattice simulations are not yet applicable 
for this regime. Our current field-theoretic knowledge in this direction comes from studying 
Nambu-Jona-Lasinio (NJL) type models with the hope that they shed light on phases of 
QCD for which the theory is strongly coupled (and gluon degrees of freedom seem to be 
irrelevant). See [2] for a review of NJL models. 

To model phenomena arising within strongly-coupled QCD, there is an emerging alter- 
native approach based on gauge-gravity duality [3]. Although it is only for a large number 
of colors Nc and large 't Hooft couphng A that one can do reliable computations using this 
duality, experience over the past decade indicates that some of these results may be applica- 
ble to QCD, which might lie within the same universality class of the theories being studied. 
One model which shares a number of common strong-coupling features with QCD, such 
as confinement and chiral symmetry breaking (^SB), is the Sakai-Sugimoto model [4], also 
known as holographic QCD. The model consists of Nc "color" D4-branes which intersect A^f 
"fiavor" D8-branes and Nf D8-branes at two (3 -|- 1) -dimensional intersections, taken to be in 
the directions. The D8- and D8-branes are located at the antipodal points of a circle in 
the direction. Imposing anti-periodic boundary conditions for fermions around the circle 
leaves the gauge bosons of the 4 — 4 strings massless but gives mass to their fermionic and 
scalar modes. At the D4-D8 and D4-D8 intersections, there are massless Weyl fermions, 
denoted V'l and V'r, which come from the Ramond-Ramond (RR) sector of the 4-8 and 4-8 
strings and transform as (Nc,Nf, 1) and (Nc, 1, Nf) of U(A'c) x U(A^f) x U(Af), respectively. 
The U(A''f) X U(A^f) gauge symmetry of the fiavor branes is the chiral symmetry for these 
Weyl fermions. At large Nc and large effective four-dimensional 't Hooft coupling Acff, in 
the probe approximation (Af -C Nc) the D8- and D8-branes are smoothly connected into a 
U-shaped configuration whose interpretation is that the dual gauge theory, which is already 
confined, is in a xSB phase. 

Moving the fiavor branes away from the antipodal points of the x^-circle leads to a 
separation in the scales of x^B and confinement. Furthermore, taking the radius of the 
circle to be infinite results in a model without confinement [5]. This non-compact version of 
the Sakai-Sugimoto model is dual to a field theory which can be referred to as the holographic 
NJL model. This model is interesting in its own right since, like the usual NJL models, it 
enables one to analyze xSB holographically in a setting where confinement has been turned 
off. 
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These models have been generahzed in a number of different ways, which have been found 
to continue to exhibit similar behavior to QCD and NJL models. For example, temperature 
can be added to the Sakai-Sugimoto model [6] and its non-compact version [7] by taking 
the D4-branes away from the extremal limit. At high enough temperatures, the preferred 
configuration is that of separated parallel D8- and D8-branes, which signals chiral symmetry 
restoration (xSR) in the dual field theory. Including a finite baryon density in the field theory 
is another generalization of these models which was studied in [8, 9, 10, 11, 12, 13]. Following 
[14], the baryons themselves correspond to S'^-wrapped D4-branes which are dissolved into 
the flavor branes and appear as their worldvolume instantons [4, 15, 16, 17]. The phase 
diagram of the Sakai-Sugimoto model at finite density and temperature has been partially 
mapped out. To some extent, this shows a resemblance to the QCD phase diagram obtained 
from analyzing field theoretic toy models. 

Despite a fair amount of success, the Sakai-Sugimoto model and its non-compact ver- 
sion have some serious shortcomings. For instance, one drawback of these models is that 
one cannot write an explicit mass term for the localized fermions because there is no direc- 
tion transverse to both the color and flavor branes along which to stretch an open string. 
Moreover, although these are models of xSB, the order parameter for spontaneous xSB is 
conspicuously absent. A promising proposal for how to modify these models in order to be 
able to compute the xSB order parameter as well as incorporate a bare fermion mass has 
been given by Aharony and Kutasov in [18] (see also [19] for related ideas). The idea is 
that, in holographic QCD and the NJL model, the left and right-handed fermions ipL and 

are localized at different points in the x-direction, so that even though ijj\^ipR is charged 
under the chiral symmetry it is not a gauge invariant operator. Hence, {ipl^i/jji) cannot serve 
as a xSB order parameter in these models. However, one can make a gauge-invariant op- 
erator out of the left and right-handed fermions by inserting an open Wilson line between 
them. The vacuum expectation value (vev) of this new operator, hereafter called open Wil- 
son line operator, or OWL operator for short, is an order parameter for xSB. It was argued 
in [18] that the vev of this operator can be obtained from the area of a Euchdean string 
worldsheet bounded by the flavor branes. In other words, this Euclidean worldsheet is the 
holographic dual of the OWL operator. The vev of the OWL operator was calculated in [18] 
for the Sakai-Sugimoto and holographic NJL models at zero temperature and density where 
a non-vanishing result was found. The proposal was then generalized in various directions in 
[20, 21], including the effects of sub-leading corrections, temperature and background electric 
and magnetic fields on the vev of the OWL operator. 

Different regions of the phase diagram of QCD are associated with different condensates. 
To determine, for example, in what regions of the phase diagram chiral symmetry is broken, 
one has to analyze the behavior of the chiral condensate (the xSB order parameter) as a 
function of temperature and density (or chemical potential) and see whether a non- vanishing 
condensate minimizes the appropriate thermodynamical potential. As we alluded to above, 
such analyses cannot easily be carried out for QCD at intermediate densities. Thus, it is 
interesting to see whether holographic models can provide some results in this direction. 

In this paper, we use the proposal of [18] to analyze how the xSB order parameter in 
the holographic NJL model responds to turning on quark (number) density, or chemical 
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potential, and temperature. Since this pertains to xSB, we expect the results obtained in 

this paper to stay qualitatively the same for holographic QCD, as well. This is in analogy to 
the case in field theory for which NJL models show, at least qualitatively, similar behavior 
to QCD in regimes where confinement is of minor relevance. 

This paper is organized as follows. In section 2, we review the holographic NJL model 
at finite temperature and quark density (or equivalently, baryon density, when quark density 
is divided by Nc). We discuss the various sources in the bulk which give rise to a density 
in the dual field theory. Along the way, we find a new configuration where S^-wrapped 
D4-branes are at a finite distance from the horizon (in fact, very close to the horizon), 
and are connected to the flavor brancs by flnitc-length fundamental strings. However, this 
configuration is unstable against small radial fluctuations in the position of the wrapped 
D4-branes, and also is less thermodynamically favored compared to other configurations 
considered in this section. 

In section 3, we calculate the xSB order parameter in the presence of finite temperature 
and density using the proposal of [18]. In the probe approximation, when it comes to 
calculating the x^B order parameter, the temperature manifests itself holographically within 
the background geometry of the nonextremal D4-branes as well as in the boundary condition 
on the Euclidean worldsheet dual to the OWL operator. The density or chemical potential, 
on the other hand, only appears through the boundary condition on the worldsheet. We 
first consider the case of zero temperature and finite quark (or baryon) chemical potential 
H. In this case, unlike the usual NJL model, the holographic NJL model cannot realize a 
xSR phase, simply because the flavor branes always remain connected. There is, however, a 
different kind of phase transition at some critical chemical potential /^cr, which is at the order 
of the baryon mass. This is a transition from a xSB phase with zero baryon density, for which 
the configuration of flavor branes is U-shaped, to a xSB phase with non- vanishing baryon 
density whose gravity dual is given in terms of a density of S^-wrapped D4-brancs attached 
to the flavor branes [11, 13]. To leading order in Acfr, these S^-wrapped D4-branes are 
pointlike [15, 16]. Working within this approximation we show that, at the phase transition, 
the xSB order parameter as a function of chemical potential is continuous, whereas its 
derivative jumps. Using a combination of analytical and numerical techniques, we show that 
for chemical potentials above but near /icr, the xSB order parameter decreases linearly with 
^ but increases for larger yU. For fi < fi^T, the order parameter is constant and equals its 
value at zero chemical potential. At finite temperature and density, the situation is more 
involved. At low enough temperature and chemical potential where the theory is in a xSB 
phase with zero density, we find that the xSB order parameter increases with temperature 
T as T^, whereas it docs not change with Increasing fi while keeping the temperature low 
causes the system to hit a critical chemical potential yUcr(T), beyond which it is in the xSB 
phase with nonzero density [11]. Our computations show that in this phase the xSB order 
parameter increases with T and initially decreases with fj, but eventually increases with both 
quantities. At high enough temperatures, the theory is in a xSR phase with non-zero density. 
By an explicit calculation involving the global structure of the background geometry, we show 
that the worldsheet dual to the OWL operator has an infinite area resulting, as expected, 
in a vanishing xSB order parameter. On the T — /j, phase diagram of the holographic NJL 
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model, there is a tri-critical point where three aforementioned phases (xSB with zero density, 
xSB with non-zero density and xSR with non-zero density) are equally thermodynamically 
favored. Combining our analyses for each separate phase, we conclude that the xSB is 
discontinuous at that point. 

2 Holographic NJL model 

We start this section by briefly reviewing the holographic NJL model and setting the notation 
that will be used throughout the paper. 

2.1 Zero temperature and density 

Consider Nc D4-branes extended in the x^^^ directions in ten-dimensional Minkowski space- 
time M^'^. To this system, add N{ D8-branes and A^f D8-branes such that they intersect the 
D4-branes at two (3 -|- l)-dimensional intersections (defects), and are separated in the x'^ 
direction by a coordinate distance (see Figure la). Focusing on the defects, there are lo- 
calized massless fermions which come from the RR sector in the 4 — 8 and 4 — 8 strings. These 
fermions are chiral with respect to the U(A''f) x U(A^f) of the flavor D8- and D8-branes, and 
transform as (Nf, 1) or (1, Nf). They are also in the fundamental of \J{Nc). The left-handed 
and right-handed fermions, ipL and ipR, interact via D4-brane gauge fields (and scalars). The 
strength of the interaction is set by the dimensionless effective coupling Aefr = As/^o- 

Stringy effects are small when gsNc is small. In this regime, when Aefr is small one can 
write down an effective action for the left and right-handed fermions by integrating out the 
gauge fields of the D4-branes, resulting in [5] 



where G{x — y, £(,) is a (4 + l)-dimensional scalar propagator. The dot in the square brackets 
in (1) denotes contraction over color indices. The above effective action indicates that the 
theory of the fermions at the defects is a NJL model where the four- Fermi interaction term 
is non-local. In the regime in which the above action is obtained, A = is the cut-off 
of the model. It was argued in [5] that, unlike the local NJL models where spontaneous 
xSB happens only above a critical coupling, chiral symmetry is spontaneously broken at 
arbitrarily weak coupling for the system of fermions described above. For large QgNc and 
large Aeff, the action (1) can no longer be trusted to provide a reliable description for the 
dynamics of the fermions. Instead, this is the regime for which gauge-gravity duality can 
be used to study the system. In order to do this, one can consider the flavor branes in the 
near-horizon geometry of Nc extremal D4-branes, described by the metric 
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(a) 




Figure 1: (a) The D4-D8 and D4-D8 intersections in R^'^. (b) At large Nc and large Aeff, the geometry of 
color D4-branes is replaced by its near horizon limit and the preferred configuration of the fiavor branes is 
U-shaped, which signals xSB in the dual theory. 



where dQ^ is the metric of the unit 4-sphere. The characteristic radius of the geometry R is 
related to the string couphng Qs and string length Ig by 

= ngsNJl = {4n)-'glNJl (3) 

where g^, is the dimensionful coupling of the low-energy theory of the D4-branes. There is 
also a dilaton and an RR 4-form flux F(4) given by 

where C(3) is a 3-form potential and e(4) is the volume form of the unit S^. The dynamics of 
the bosonic degrees of freedom of the fiavor D8-branes are determined by the Dirac-Born- 
infeld (DBl) plus the Chern-Cimons (CS) action 

S = 5dbi + ^cs = -1^8 J d^C e-^TT^-det{g + 27ra'F) + ^ ^ C A Tr e^-"'^. (5) 

Fab s-nd gab in (5) are the field strength and the puUback metric, respectively. The D8-branes 
(as well as the D8-branes) are extended in the {t, x, S^} directions and form a curve u — u{x'^) 
such that the asymptotic distance between the D8- and D8-branes is Eq: ?i(±£o/2) = ma. 
Since the dilaton grows with u, ma cannot take arbitrarily large values. Indeed, in order for 
the supergravity approximation to be valid, u\ has to satisfy ua <^ a'N^I^/gl [22]. Since 
the D8-branes are in the background of Nc D4-branes where only C(3) is turned on, the 
Chern-Simons part of the action reads 



'S'cs 
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where uj^5){A) = TryAF'^ - jqA^j is the Chern-Simons 5-form satisfying duj^5){A) = 

Tr F A F A F, and the integration on the right hand of (6) is over a five-dimensional space 
A4^ X R+ spanned by the and u directions. 

Analyzing the equations of motion coming from (5), one can set the gauge fields equal 
to zero. The equation of motion for u = u{x^) then shows a U-shapcd profile for the 
fiavor branes [5]. This U-shapcd solution, which is energetically favored, is interpreted as 
representing xSB, where at large u the U(iVf) x U(A^f) symmetry is manifest while at the 
lowest position of the U-shapcd profile there is just one U(A^f) factor (see Figure lb). Thus, 
the holographic NJL model at large Nc and at strong coupling, Aefr ^ 1, has a vacuum in 
which chiral symmetry is broken. 



2.2 Finite temperature and density 

At finite temperature, one considers the fiavor branes in the near horizon geometry of Nc 
non-extremal D4-branes whose metric, in Lorentzian signature, is 

where Ut is the radius of the event horizon, T is the temperature of the black D4-branes, 
and R is given by (3). The relationship between T and ut in (8) can easily be understood 
by Wick rotating the metric in (7) to the Euclidean signature and demanding the geometry 
to be regular at ut- The dilaton and the RR 4-form fiux are the same as for extremal 
background, and are given by (4). 

At finite temperature, the analysis of the DBl action for the fiavor branes shows two dis- 
tinct profiles: parallel and U-shaped. The parallel profile represents fiavor branes separated 
from each other going down to the horizon of the background geometry. The field theory 
interpretation of this solution is that the U(A^f) x U(A^f) chiral symmetry of the holographic 
NJL model is unbroken. The U-shaped profile, on the other hand, represents connected fia- 
vor branes which stay above the horizon. The dual interpretation of this solution is that the 
holographic NJL model is in a phase with spontaneous xSB. There is a critical temperature 
of T-^^sB — 0.15 (in units of £q^) below which the U-shaped solution is the thermodynamically 
preferred one^. For temperatures above T^sb, the parallel solution is favored, which results in 
the restoration of the chiral symmetry of the holographic NJL model. This phase transition 
is first order [7]. 

The U(1)b baryon symmetry in the holographic NJL model is associated with the U(l)v' 
subgroup of the U(iVf) x U(iVf) ~ SU(iVf) x SU(iVf) x U(l)y x U(1)a symmetry of the 

^Strictly speaking, there are two U-shaped solutions at small enough temperatures. One solution bends 
closer to the horizon while the other one stays farther away from it. It turns out that the solution which 
stays farther away from the horizon is always thermodynamically favored [7, 11, 23]. 
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flavor branes. To be more speciflc, the baryon number is deflned to be A^~^ times the 
charge of U(l)y C U(iVf)v'. The timehke component of the U(l) gauge field on the D8- 
and D8-branes couples asymptotically to ip^j^ipL and ip\^tpR, respectively. Thus, introducing 
a chemical potential /i for the baryon number operator Qb = N~^{ip]^ipL + '4'r4'r) i^i the 
holographic NJL model amounts to turning on the timelike component Ao{u) of the U(l) 
gauge field on the fiavor branes and requiring that Ao(oo) — /i on both branches. More 
precisely, according to the standard gauge-gravity dictionary, the non-normalizable mode of 
Ao gives the chemical potential while the normalizable mode is interpreted as the baryon 
(quark) number density. In the rest of this paper, unless there is confusion, we will use the 
baryon number density ub and the quark number density riq interchangeably, and sometimes 
refer to them simply by the density. It is understood that ub — N~^n^. 

It is convenient to work with dimensionless quantities. We make {t, a;, u} as well as Aq{u) 
dimensionless by scaling out factors of R and 2T:a' according to 

t X _ u 2'Ka' 

-^i, R^'^' -^Aq{u) ^ Ao{u). (9) 

With the above rescahng, the metric takes the form 

^ = G^^dx^dx'' = ( - f{u) de + dx^) + u- 'i (-^ + u^dni) , (10) 

and T is now the dimensionless temperature given by 

(11) 

Also, to make the notation less cluttered, we denote the x^ direction by Due to symmetry 
considerations, we choose a purely radial ansatz for the timelike component of the U(l) 
gauge field: Aq — Aq{u). With fiavor branes forming a curve u{x) in the x — plane and 
with ^0 turned on, the action for the total configuration takes the form 

S = -2N,C r duJ ^/(«)«V + [!-<], (12) 

Jut 

where the primes on x Aq denote derivatives with respect to u, and Uf is the lowest 
radius to which the fiavor branes descend. Also, the dimensionless constant C is given by 

C = ^i?Vol(R^)Vol(S^), (13) 
9s 

where Vol(R^) is the dimensionless volume of space spanned by the {t ^} directions. 

Varying (12) with respect to Aq and x results in the following first integrals of motion: 

d = / (14) 
^Jf{u)u^x" + [1 - ^0 ] 

c. = , , (15) 

^f{u)u^x" + [1 - ^ ] 
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where d, in some units, is equal to the electric displacement on the the flavor branes and 

is a parameter which determines the shape of the flavor branes in the x ~ plane. It 
is sometimes convenient to trade A'q in the action (12) with d by performing a Legendre 
transformation on the action. In so doing, we obtain 

S^S- 2NfC j du A'^d = -NiC j du J J (l + f{u)u^x'^) (l + u'^dA , (16) 

Juf Jut y 



where we made use of (14). Solving for in (14) and expanding the result in the asymptotic 
region u ^ oo yields 

^(«) = Ao(oo)-^^ + ---, (17) 

which shows that d appears as the coefficient of the normalizable mode oi Aq. d is propor- 
tional to the baryon number density tzb = {Qb) (of) the quark number density riq) while 
Aq{oo) equals the dimensionless chemical potential /i for the density. 

We start the analysis of the above first integrals of motion (14) and (15) by first consid- 
ering the parallel embedding of the flavor branes for which x' = 0, or equivalently = 0. 
In order for the induced metric on the flavor branes not to change sign in this embedding, 
^0 should satisfy an upper bound of (^q)^ < 1. The physical meaning of this bound is 
straightforward to understand: since F^t — Aq is the radial electric field on the flavor branes, 
the bound simply states that Fut must be less than a critical value in order to prevent 
open strings on the branes from breaking up into pairs. Setting = in (14), there is a 
trivial solution Aq{u) = constant for which d = 0. Since the near-horizon geometry of Nc 
non-extremal D4- branes given by the metric (10) is static, the surface u — ut is a, Kilhng 
horizon. Indeed, it is a bifurcate Killing horizon at which the timelike Killing vector dt 
vanishes [24]. This, in turn, implies that must vanish at the horizon: Aq{ut) = 0. Thus, 
this conflguration of flavor branes represents a chirally-symmetric phase of the holographic 
NJL model with no density turned on. In order to obtain a nontrivial density, one must turn 
on a nontrivial (radial) electric field F^t = A'q whose fines end on a density of sources. 
For a parallel embedding, the fiavor branes intersect the horizon, and hence the electric field 
lines on the branes can end there. Thus, one can think of the horizon as providing a density 
of sources causing a nontrivial profile for the radial electric field on the branes. Solving (14) 
for Aq for the parallel embedding, we obtain 

p d d 
M'^) = / / , , du, = / du, (18) 

Jut VU + " Jut V 1i + « 

where we set Ao{ut) — 0. Note that as given by (18) satisfies {A'^Y < 1 for it > ut- If we 
work in the grand canonical ensemble where T and fi are kept fixed, we then need to compare 
the grand free energy Q{T,fi) of the two above-mentioned parallel solutions (one with zero 
density and the other with non-zero density) to see which one is thermodynamically favored. 
Holographically, 0(T, /x) is given by the Euclidean on-shell action of the configuration of 
the fiavor branes, appropriately renormalized. For parallel embeddings, ^{T, /i) equals the 
Euclidean continuation of (12), once x' = and (18) is substituted in (12). Since both 
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of these parallel embeddings have the same UV divergences, one can just subtract their 
Euclidean on-shell actions, without bothering to renormalize them, and get a finite result. 
Using (18), one obtains that for all T and /j,, the parallel configuration with a non- vanishing 
density is always favored. 

For an embedding where the branes connect, one has x' 7^ 0. There is a solution for 
which x'(m), for half of the configuration, say, where x ^ [0,^o/2], is a smooth single-valued 
function of u. In this case, %'(«) is given by (15), and ^o(^) ~ 0- This solution represents a 
xSB phase of the holographic NJL model with zero density. To be more precise, for io less 
than a critical value there are two such solutions, corresponding to two distinct values of c^. 
To study the holographic NJL model at finite density, an electric field on the flavor branes 
must be turned on. To better analyze this case we use the u = u{x) embedding for which 
the action for the flavor branes reads 

reo/2 I 

S = -2N,C / dxu-^JfiuW + [{dx^r - [d^Aof]. (19) 
Jo 

The flrst integrals of the equations of motion are 

d = "^^^^^"^ (20) 



Cx = / ==• (21) 



As we alluded to earlier, in order to turn on a baryon chemical potential in the holographic 
NJL model, the asymptotic value of Aq must be the same on both branches of the flavor 
branes [13] which, by virtue of symmetry, means that Aq is symmetric about x = 0. Assuming 
an everywhere smooth d^Ao, one deduces that at x = is an extremum, i.e. d^Ao\^=o = 0. 
From (20), this yields d = with Aq being a pure gauge everywhere on the branes. In order 
to avoid this trivial outcome, one can assume that O^Aq is discontinuous at x = [12], which 
means that there is a density of charges located at x = providing a source (or a sink) for 
the radial electric field lines. In the presence of this source, a smooth embedding of the fiavor 
branes can no longer be at equilibrium simply because there is no tension in the M-direction to 
counteract the gravitational force exerted by the source. Thus, for a non-vanishing density 
we discard smooth embeddings as they do not represent stable configurations. A stable 
configuration may be obtained if the connected embedding has a cusp at some radial point, 
say at M = ut, in which case Ut and the opening angle of the cusp can be determined by 
balancing the forces at the cusp. 



Point-like D4-branes One possible source is fundamental strings attached, at one end, 
to the fiavor branes at u = Ut- To obtain a net charge, the string configuration has to 
connect at the other end to something other than the flavor branes. One possibility is that 
the strings end on D4-branes to form a baryon vertex in the background geometry. These 
D4-branes are extended in the {i, S^} directions and are at flxed points in the rest of the 
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directions. The location of the D4-branes in the u direction can be obtained by extremizing 
the on-shell action. Equivalently, if we assume that the D4-branes are at some u = Uc, then 
Uc can be determined by requiring that the total tension in the u direction at Uc vanishes^. 
We denote the number of these D4-branes by 2iVD4 and their density by 2nD4 (where the 
factors of 2 are just for later convenience) such that A^d4 = nD4Vol(R^). Since the gauge 
fields on the worldvolume of these D4-branes play no role in our discussion, it is consistent 
with the equations of motion (obtained from the DBI action) to set them equal to zero. The 
DBI action of the D4-branes, localized at Uc {ut < Uc < Ut), then reads 

So. = -A. Jd^ie-^ Tt y^d^ = -2A-o.no. / du S(u - u,)u/f(u), (22) 
where gj^^ is the induced metric on the D4-branes and K^^ is a dimensionless constant given 

by 

Xd4 = — i?'Vol(R^)Vol(S^). (23) 



9. 
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There are Nc strings attached to each D4-brane giving a total of 2A^f = 2NcNd4 strings 
with a density of 2nF = 2A^c?^d4- Again, the factors of 2 are for convenience in subsequent 
calculations. The action for 2A^f strings coupled to Aq{u) on the flavor branes at Ut and 



attached to the D4-branes at Uc is 



Sp = J dV^-detc/F + 2iVF j A^dx^ 

= -2KpNcnYii{ut-Uc) + 2KYNcnjii J Ao{u)5{u - Ut)du, (24) 

where (/p is the induced metric on the worldsheet of the strings and Kp is a dimensionless 
constant given by 

i^F=^Vol(R^). (25) 

Note that A^ and x^ in the first hne of (24) are dimensionful, whereas in the second line their 
dimensions have been scaled out according to (9). We denote the tension of the D4-branes 

in the u direction at -u = -u^ by /d4|«c- '^^^ compute it by first varying the on-shell 

action of the D4-branes with respect to Uc and then multiplying the result by a factor of 
[Guu{uc)]~^^'^ to account for the proper tension [11]. The tension of the strings at m = Uc, 
denoted by /f|uc) can similarly be computed. The magnitude of these tensions takes the 
form 

3 

/d4|«, = KD4nB4.Uc[S - f{Uc)], (26) 



/fU, = 2i^FiVenD4«ly7K), (27) 



^This statement assumes that the D4-branes, as well as the strings, do not interact with one another, 
and that the strings are all stretched in the u direction. 
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where to obtain (26) and (27) we have assumed that f{uc) is non-vanishing. Using (3) along 
with Vol(S^) = |7r^ = 27r(27rZs)-^ and h = Va', one deduces that 3^04 = NcKp. 

It is easy to see that, for ut < Uc < Ut, the two tensions given in (26) and (27) are 
equal only at Uc — 1.08-ut- One can show that, when perturbed away from Uc — I.OSmt- 
the D4-brancs will accelerate further away. This indicates that the D4-brancs arc not at a 
stable equihbrium at Uc — I.OSmt- A stable configuration can potentially be obtained if the 
D4-branes are at Uc — Ut- For Uc — Uf, the strings have zero tension, whereas the tension 
of the D4-branes is non- vanishing. Thus, in order to obtain a stable configuration, the non- 
zero tension of the D4-branes must be equal to the radial tension of the fiavor branes at 
Ut- Another potentially stable configuration is obtained when the D4-branes cross through 
the horizon, so that a bundle of strings stretch between the fiavor branes and the horizon. 
Again, for this configuration to be stable, the tension of the fiavor branes at Uf must equal 
the tension of the strings. 

The action for the total configuration, which includes 2A^f fiavor branes and a density of 
2nD4 D4-branes located at Ut, reads 



S = -2NiC duulJf{u)u^x'^ + [l-A'^^^]-2K^inB4 duS{u-Ut)u^/J^ 



Setting the variation of the above expression with respect to Aq equal to zero, we obtain 



In order for the configuration to be stable, the tension of the 2A^d4 D4-branes at Ut must be 
equal to the tension of the flavor branes. We denote the tension of the flavor branes in the u 
direction by /flavor |u- Due to the symmetry of the configuration along the w-axis, the forces 
in the x direction are already balanced. Keeping £0 fixed, the condition for the forces to be 
equal at Ut is 



where is given by (15). Note that ^04 in (30) is related to d through (29), and is related 
to io. Thus, for fixed io (which is our assumption), (30) relates Ut to d. 

Fundamental strings As stated before, the configuration of the fiavor branes with a 
density of 2nF = 2NcnY)4 strings stretched from the fiavor branes ai u = Ut to the horizon is 
another potentially stable configuration which can source d. The action for this configuration 
is 




(28) 



d = {NiCy^KpN^riBA- 



(29) 



2N,Cu;^^^f{ut){ul + uld^) - 4 = nD4i^D4«|[3 - /(«*)], 



(30) 



S ^ - 





(31) 
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where the orientation of the strings is chosen to be upward, i.e. from the horizon to the 

flavor branes. The relationship between the electric displacement d and rip is the same as 
the one in equation (29): d = {NfCy'^KpriF = {NfCy'^KpNcriBA- 

Keeping £o fixed, the condition for the forces to be equal at Ut reads 

NiCuf' ^Jf{ut){ul + uld^)-cl = KpUFuj V/K). (32) 

Note that, due to the symmetry of the configuration along it-axis, the tensions in the x~ 
direction are already balanced. Given (29) and the fact that is related to £o, (32) relates 
Ut to d for the configuration of the fiavor branes and the fundamental strings. 



Instantons (dissolved D4-branes) Another possible source for d is a density 7204 of 
D4-branes^ which are dissolved into the D8-branes and appear as worldvolume gauge fields 
carrying instanton charge 



nD4 = ^ / TVF A F. (33) 



The integral in (33) is over a four-dimensional space A4 spanned by {x^,x'^,x^,u}. These 
instantons couple to the Abelian gauge field on the worldvolume of the D8-branes through 
the Chern-Simons part of the action. The electrostatic repulsion among the instantons tends 
to spread them out and give rise to a distribution of instantons in the w-direction, whereas 
the gravitational force acts in favor of localization and makes the instantons point-like. 
The competition between these two forces gives the instantons roughly an effective size of 
p ~ ^eff^"^ [^^' 16]. Note that, for A^f = 1, finite size instantons cannot be topologically stable 
and tend to become point-like, whereas for A^f > 1 this is not the case. Regardless of the 
value of Af, the instantons can effectively be treated as point-like objects for large enough 
Aeflf. Assuming a uniform density u-da of D4-branes localized at u — Ut, the Chern-Simons 
term (6) reads 

Scs = ;^nD4A^cVol(R^) / Ao{u)5{u - ut)du, (34) 

which, in turn, results in ci = A'fA'c^d4- Note that 77,04 is the density of D4-branes needed to 
source d for half of the fiavor configuration. The actual density of instantons for the entire 
configuration of the D8-D8-branes should therefore be 2nD4. In the point-like instanton 
approximation, one can show that the DBI plus Chern-Simons action for the configuration of 
the flavor branes with a density of (point-like) instantons on their worldvolume is equivalent 



^Since D6-braiies can also end on D8-branes, a bundle of D6-branes stretched between the flavor branes 
and the horizon is another stable configuration providing a source for d. Since a D6-brane ending on D8- 
branes is a magnetic monopole on the worldvolume of the D8-branes, such a configuration will source a 
density of magnetic charges in the boundary theory. We do not consider such configurations here because 
we are mainly interested in turning on sources for objects which are electrically charged under Aq, such as 
quarks. 
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(up to an term) to the action (28) for the flavor branes with a density of point-hke D4- 

branes locahzed at Ut [15]. We work in this approximation in this paper and can therefore 
use these two configurations interchangeably. Going beyond this approximation, it turns 
out that the configuration of fiavor branes with a density of finite-size instantons on their 
worldvolume is thermodynamicaUy favored over the configuration of flavor branes with a 
density of point-hke D4-branes [13]. 



3 xSB order parameter 

In the holographic NJL model, the left and right-handed fermions iJjl and iJjr are localized at 
different points in the x direction^. At strong coupling, X^g ^ 1, the model is essentially flve- 
dimensional. Thus, the one-point function of the chiral operator '0l V'iR cannot be identifled 
as a xSB order parameter simply because it is not gauge invariant. Here, ■ ■ ■ denote 
fiavor indices. One can make a gauge-invariant chiral operator by inserting an open Wilson 
line between ipL and ipR. This leads to the so-called open Wilson line operator [18] 



Oiixf^) = i^l{x\x= -WV exp / {lA^ + ^)dx ^iR(a;^x = £o/2), (35) 



where A-^ is the component of the gauge field in the x direction. The non-local operator 
0] (x'*) depends on the choice of the contour for the open Wilson line insertion. For different 
contours, one obtains different operators of the above type. The simplest operator is the 
one for which the contour consists of a straight line traversing the x direction. The operator 
0'^j{x^) indeed represents a class of non-local chiral operators whose one-point functions can 
potentially serve as xSB order parameters of the holographic NJL model (or holographic 
QCD). One can slightly generalize 0](x'*) by considering ■0^, and '4>r to be located at two 
different points in the x'^-directions, in which case we denote the corresponding operator 
by Oj{xi,X2)- In analogy with the well-stablished holographic dictionary for Wilson loops 
[26, 27], it was conjectured [18] that Oj{x'^), or in general C*(a;^, Xg), is holographically dual 
to a Euclidean string worldsheet bounded by the fiavor branes as well as the contour of the 
open Wilson line Bit u — ua, while being at fixed points in the rest of the directions^. More 
precisely, the one-point function of the operator is given by 

{O]) = 5]e-^-, (36) 

where is the regularized action of the above-mentioned Euclidean worldsheet. This world- 
sheet has Neumann boundary condition where it attaches to the fiavor branes and Dirichlet 
boundary condition on the contour aX u = u\. Unless there is potential for confusion, we 



"^This is also the case for the Sakai-Sugimoto model and, in general, for transversely-intersecting D-brane 
models considered in [25, 23]. 

^Thc authors of [19] discuss siinihu- ideas on modifying the holographic QCD and NJL models in order to 
be able to introduce a mass term for the fermions and compute chiral condensates. Refs. [28, 29] (see also 
[30]) take another approach for computing fermion masses and chiral condensates in these models, which is 
based on including the dynamics of open string tachyons in the DBI action of the flavor branes. 
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will delete the x'^-dependence of Oj{x^) as we already did in writing equation (36). Also, 

we sometimes omit the flavor indices on the operator and write it simply as O. Figure 2a 
shows the worldsheet dual to the simplest operator of the type in (35), the one for which 
the contour C is a straight line extended in the x direction at u = u\. Figure 2b depicts the 
worldsheet dual to an operator with a curved contour C. In practice, it is hard to find the 
one-point function of the Oj{x^) operators dual to such curved worldsheets. Calculating the 
one-point function of the simplest open Wilson line operator using the conjecture (36), one 
then obtains a non- vanishing result for (Oj) in the vacuum (a xSB phase) of the holographic 
NJL model (or holographic QCD) at zero temperature and density [18]. In this paper, we 




Figure 2: (a) The worldsheet dual to the simplest open Wilson line operator: the one for which the contour 
C at u = Ma traverses a straight line in the x direction and has a length of Iq. (b) A worldsheet whose dual 
open Wilson line operator has a curved contour C. 

focus on the aforementioned simplest open Wilson line operator and calculate its one-point 
function, the xSB order parameter, in different phases of the holographic NJL model at 
finite density and finite temperature. We would like to know how the xSB order parameter 
responds to simultaneously turning on chemical potential for the baryon number density and 
temperature in the holographic NJL model. We start our analysis by considering first the 
case of zero temperature but finite density of the model. 

3.1 Finite density and zero temperature 

At zero temperature but finite density, the configuration of fiavor branes and fundamental 
strings should be discarded. This is because the background geometry does not have an 
event horizon, in which case the strings have to be stretched from the flavor branes at 
u — Ut all the way down to it = 0, where the background curvature (in units of Is) becomes 
large and the supergravity approximation breaks down. Thus, we are left with either a 
U-shaped conflguration of flavor branes for which = 0, or else a conflguration with flnite 
density for which point-like D4-branes are attached to the flavor branes at Ut and source 
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the finite value of d. On physical grounds, for small enough values of chemical potential, 

one expects the baryon states not to be able to populate, which results in the solution 
with zero density being preferred. Beyond a critical chemical potential, which is basically 
equal to the mass of the baryons (in the approximation where the interactions amongst the 
baryons are neglected), the expectation is that the solution with a non- vanishing density 
is thermodynamically favored. It has been shown in [11, 13] that such expectations are 
met in the Sakai-Sugimoto model: for fj, < /Xcr, where in the holographic setup /icr is the 
total mass of the point-like D4-branes, the U-shaped flavor configuration has lower free 
energy, whereas for /j, > the configuration of fiavor branes with a density of point-like 
D4-branes localized at Ut is favored. It is easy to show that the holographic NJL model 
also shows similar behavior. Although the holographic NJL model at zero temperature and 
finite density behaves qualitatively the same, the field-theoretic NJL models can exhibit 
rather different behavior. In particular, at zero temperature, the field-theoretic NJL models 
undergo a xSR phase transition at high enough baryon chemical potentials, or densities, (for 
example, see [2]). 

On the other hand, this is not what occurs in the holographic NJL model up to leading 
order in large N^., since the flavor branes are always connected at zero temperature. Recall 
that the actions of both point-like D4-branes and flavor branes scale as Nc (to be more 
precise, NfNc), whereas the action for the background geometry scales as A*"^. While the 
configuration of fiavor branes and point-like D4-branes does not back-react on the geometry 
to leading order up in N^, it does so through the corrections. Taking into account the 
back-reaction, it is plausible that the holographic NJL model may exhibit a xSR phase at 
zero temperature but finite chemical potential. It is interesting to study this particular issue 
further. In this section, we calculate the xSB order parameter (0*) up to leading order in 
large Nc and large \cs for the above two phases of the holographic NJL model and determine 
its dependence on both the baryon chemical potential fi and baryon number density ~ d. 

In order to calculate (0]), we need to know the profile u = u{x) of the fiavor branes. 
Recall that the worldsheet is bounded by the fiavor branes defined by the curve u = u{x) 
and by the contour at u — ua- For n < /Xcr, the configuration of the flavor branes is U-shaped 
where the density is zero. Setting f{u) — 1 and dy^Ao = in (19), one obtains from the 
action the flrst integral of motion 



where uq is a non-vanishing constant for a U-shaped conflguration. Given that u{±io/2) — 
u\, (37) yields 



where to obtain (38) we have ignored terms subleading in u\. We want to calculate the area 
of a Euclidean worldsheet bounded by the curve u = u{x) which is given by (37). In the 
absence of gauge flelds on the world volume of the flavor branes, the string action is given by 




(37) 




(38) 





(39) 
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The string equation of motion then reads 



dadpx^ - Tlpd^x^ + V^,^dax''dpx'' = 0, (40) 

where 7 = {(7^,0"^} denote the worldsheet coordinates while are spacetime coordi- 
nates, and are the connections that correspond to the induced worldsheet metric 
ha[) and the background metric G^,^, respectively. It is straightforward to see that taking 

(7^ = Rx, = Ru, and x'^{a^, a^) = {[i^ x, u), (41) 

where Xq is a constant, solves the string equation of motion (40). This worldsheet, depicted 
in Figure 2a, is the simplest worldsheet and was considered in [18]. The on-shell action of 
this worldsheet takes the form 

?2 



Sw = / dudx\/ GuuGx^x'^ 



R 



27ra 



where U\ is a cutoff. The linearly diverging piece in the third line of (42) has no effect on 
our discussion and can, in fact, be absorbed in the definition of (Oj) [18]. This term could 
either be subtracted away by hand or else done so automatically by performing a Legendre 
transform on the string action along the lines of [31] (see, also [32]). Dropping the linear 
divergent term and using (3), (38) and the fact that A5 = {2T^YgsNcls, we arrive at 

= -cAeflf, (43) 

where c = ^b(^^, \^ /128 ~ 0.008. Note that (43) is the area of the worldsheet up 

to leading order in X^q. The subleading corrections in Aes for this worldsheet have been 
calculated in [20]. In this paper, we determine (C*) up to the leading order in Aeg. For 
the xSB order parameter in a phase of the holographic NJL model where chiral symmetry 
is broken and the density is zero (or equivalently, for fi < /icr), substituting (43) into (36) 
yields 

= Sy^^«. (44) 

Not surprisingly, the result in (44) is the same as the one obtained in [18] for the holographic 
NJL model at zero temperature and zero baryon chemical potential. 

For d 7^ 0, there is a density of point-hke D4-branes on the gravity side attached to the 
configuration of fiavor branes at Ut- To find the profile of the fiavor branes in this case, we 
eliminate d-^Ag between (20) and (21) which, having set f{u) — 1, results in 

{d^ur = -fu' + u'd'-ul), (45) 
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where ul = is a constant. Prom (45) we obtain 



1^0 4 

— = 



and from (30) we have 



du 



(46) 



(47) 



We are not aware of a closed-form solution for the integral in (46). However, one can perform 
an analytical approximation in the limit of small d. Namely, the integrand in (46) can be 
expressed as an expansion for d^ /u^^ <^ 1. Defining = + u^d'^, (47) can be approximated 
to leading orders to be 
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where 



In{e) 



dv 



l+e2 V 



and 



72 u^' 



(48) 



(49) 



Note that, to obtain (48), we have dropped terms subleading in u\. Expanding the integrals 
in (48) and keeping terms up to yields 
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(50) 



For fixed £o, (50) along with (47) provides a relation between Uf and d, for small values of 
d. As is well known, in holographic models the radial direction translates into an energy 
scale of the dual field theories. In particular, in our holographic setup, Ut can be thought 
of as the energy scale at which spontaneous xSB happens. It can be shown from (50) that 
this scale decreases with the density d in the limit of small d. Numerically, one can verify 
that this scale keeps decreasing with the density d until ~ 0.05, beyond which it increases. 
Indeed, this behavior has already been observed in [11], where a plot of Uf versus density 
was presented for the Sakai-Sugimoto model at finite temperature and density. 



For 0? 7^ 0, the string action (39) now has an additional boundary term 



b.d. 



(51) 



which modifies the string equation of motion (40). But, it can easily be shown that the 
worldsheet extended in the x — plane and at a fixed point in the rest of the directions 
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is still a solution. Having subtracted the linear divergent term and ignoring terms that are 
subleading in ux, the area of this worldsheet reads 



K 4 



J * 



du 



u 



(52) 



Changing variable to 
mated by 

S ~ 



u 



u d and working in the d /ul -C 1 regime, (52) is approxi- 
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Ehminating u* between (48) and (53) yields, to leading order. 



This results in 



5w = -cAeff(l-2.1£^rf 



{0))d = m 1 - 0.02Aeff£^rf 



(53) 



(54) 



(55) 



where denotes the xSB order parameter at finite density. In order to obtain 

for arbitrary values of rf, the integrals in (48) and (53) should be evaluated numerically. 
Figure 3 shows how (normalized by its zero density counterpart (Cj)) varies with 

the density d. Although not shown here, it could easily be verified that the scale of xSB, 
which is proportional to Ut, exhibits rather similar behavior as a function of the density d 
as what is shown in the plot of (C'})ci in Figure 3. Namely, they both decrease with d until 
d ~ 0.05, beyond which they both increase. It has previously been found in [21] that the scale 
of xSB and the x^B order parameter {Oj)d also share similar behavior in the presence of 
background electric and magnetic fields in the holographic NJL model: they both increase 
(decrease) with a magnetic (electric) field. However, these quantities only have the same 
qualitative features at zero temperature. In particular, in the holographic NJL model at 
finite temperature, ut — ut decreases with temperature, whereas the xSB order parameter 
increases with temperature. 

Given the relationship between the density d and the chemical potential fi, we can use 
(55) to determine the leading-order dependence of the xSB order parameter on the chemical 
potential. It was shown in [11, 13] that, in order to identify Aq{oo) with the chemical 
potential, one must identify Aglut) with the total mass of point-like D4-branes localized at 
Ut (the baryons). Rescaling the chemical potential according to {2NfC)ii — > /i, then from 
(14), (15) and (22) one obtains 



U 



du + -Ut, 



■y^u^ + u^d^ - ul 37rQ; 



(56) 
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Figure 3: {0*j)d/{0^j) as a function of d. We have set R = Iq = Svra' = 1, which all together gives 
Aeff = Stt^. We have also chosen the cutoff to be ma = 10. The plot is qualitatively the same for other values 
of io and ua- Although we cut the plot at d ~ 0.3, is monotonic in d for d > 0.3. 



where Uf is given in (47). The value of for d = 0, denoted by fi^r, is 
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(57) 



For n < Her, the U-shaped configuration of the fiavor branes is dominant, whereas for /x > Her 
the configuration with lower (grand) free energy is made of fiavor branes with a density of 
point-hke D4- branes localized on the fiavor branes at Uf. For (P/ul <S 1, the expression in 
(56) approximately equals 

2v^i?2 
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45 ^ , , 2 
e+y^7(e)e' 
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16' 2 



Qttq; 

where in the second line we have kept terms up to order 
(58) can be rewritten as 

Substituting (59) into (55) gives 



— e , (58) 

dna' / J' ^ ^ 

Using (48), the expression in 



(59) 



(60) 



for < — Her -C Her- Notc that (0})// denotes the xSB order parameter as a function of 
chemical potential. Thus, (O*)^ decreases with fi for small enough /i — ficr (compared to 
Her)- For arbitrary /i, {^j)iJ. can be calculated numerically. The result is given in Figure 4, 
where we have again taken R — — 27ra' — 1 and u\ = 10. For this choice of parameters. 
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~ 0.35. Although Figure 4 is for particular values of the parameters, 

we find that the qualitative behavior of {Oj)^ is the same for other values of io and u\. Note 
that at the phase transition which happens at /icr — 0.35, the xSB order parameter as a 
function of chemical potential is continuous, whereas its derivative jumps. 
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Fi gure 4: {Oj)^/{Oj) as a function of /x. For /x < Hcr — 0.35, the U-shaped configuration is dominant and 
{Oj)n does not change with fi. For fj, > ^ct, the configuration of fiavor branes with a density of point-Uke 
D4-branes is dominant and (C])/* initially decreases linearly with for — /icr Mcr- 

3.2 Finite density and temperature: xSB phase 

As we alluded to earlier, at finite temperature but zero density there is a first-order xSR 
phase transition at the temperature T^sb — 0.15 (in units of -^q Below T^sb, the U-shaped 
solution is energetically favorable and the interpretation is that the chiral symmetry of the 
holographic NJL model is spontaneously broken. For temperatures greater than the critical 
temperature, the configuration of parallel branes is thermodynamically favored, which is to 
say that chiral symmetry gets restored. It was shown in [21] that the xSB order parameter 
increases with temperature as the critical temperature is approached from below. In the 
xSR phase, on the other hand, the xSB order parameter vanishes since the dual worldsheet 
has infinite area. 

At finite temperature and density, the critical temperature above which chiral symmetry 
is restored decreases with the chemical potential [9, 11]. At small enough values of tem- 
perature and chemical potential (compared to and /icr(T), respectively), the U-shaped 
configuration of fiavor branes is favored so that the dual field theory is in a xSB phase with 
no density. Increasing the temperature (while keeping the chemical potential fixed at a low 
value), one finds that parallel branes supporting a non-trivial density are favored. This im- 
plies that the dual field theory is in a xSR phase at finite density. Increasing the chemical 
potential (and keeping the temperature fixed at a low value) results in fiavor branes with a 
density of point-like D4-brancs being favored. Thus, the model is in a xSB phase at finite 
density. At high enough temperatures and chemical potentials, the configuration of parallel 
flavor branes which supports a non-trivial density is preferred, implying that the dual field 
theory is in a xSR phase at finite density. A schematic diagram of the various phases of 
this model with respect to the temperature and chemical potential is given in Figure 5. A 
similar diagram has already been obtained numerically in [11] for the Sakai-Sugimoto model. 
The only difference here is that, for the holographic NJL model, there is no confinement- 
deconfinement phase transition. 
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Figure 5: A schematic phase diagram of the holographic NJL model at finite temperature and chemical 
potential. 



We will denote the xSB order parameter at finite density and finite temperature by 
(C'})d,T- In order to calculate (O*)^,^ at finite temperature and density, we note that the 
profile of flavor branes satisfles 



1 + u-^d^ 
l + u-%d^u)^ 



(61) 



where is a constant. Given the boundary condition u{±£q/2) = ua, one obtains from (61) 
(ignoring terms subleading in ua) that 



4 



du 



2 "*Ju, uy^^f{u)[f{u){u^ + uW)-uiy 
where Ut (which satisfies (30)) is implicitly given by 



(62) 



ul = f{ut)[ul + u',d'] - -utd'[3 - fim)]'. 



(63) 



It is easy to show that the string equation of motion (40) has a solution with a worldsheet 
extended only in the x — w plane and at a fixed point in all other directions. The regularized 
area of this worldsheet (with the linear divergent term subtracted) reads 



'S'w — 



du 



uy^^/f{u)[f{u){u^ + u^d^) - ul 



(64) 



At low temperatures and small densities, the configuration of flavor branes with a density 
of point-like D4-branes is dominant. This results in the holographic NJL model being in 
a xSB phase at finite density. In this regime, d'^/ul <^ 1 and ut/u^, <^ 1. Defining — 
f{u){u^ + u^d^) enables one to approximate io and in (62) and (64) as 
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(65) 
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and 
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T 



respectively. Eliminating between (65) and (66) gives 

= -cAeff [1 - 2.Ul d+ 136.2£^ T^] , 



(66) 



(67) 



where we have used (11) to express ut in terms of the dimensionless temperature T. This 
results in 



(0})d,T = {O]) [1 - 0.02Aeff£^ d + 1.09Aeff£^ T^] . 



(68) 



Thus, for small values of temperature and density, the xSB order parameter of the holo- 
graphic NJL model decreases (linearly) with density and increases with temperature as T^. 

The relationship between /i (rescaled by 2NfCii — > /i) and d at finite temperature reads 



u 



zdu + 



Ut 



^f{u){u^ + u^d^)-uf^ ' Sttq;'' 
where Ut is given in (63). The value of // at o? = and low temperature is 

-7 



(69) 
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(70) 



from which we see that, at least in the limit of small temperatures, the critical value of fi 
increases when the temperature is turned on. For fi < ficr{T), the U-shaped configuration 
is dominant, whereas for fi > iXcriT) the dominant configuration is made of fiavor brancs 
with a density of point-like D4-branes. Figure 6 is a numerical result which shows how 
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Figure 6: Plot of as a function of T. 
ljLcr{T) depends on the temperature. Although Figure 6 shows HcriT) dramatically decrease 
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for larger temperatures, this does not actually occur since this would only take place for 
temperatures which are greater than where the x^R phase transition happens. 

Since a non-zero temperature does not alter the relation (59) between fi and d at the 
lowest order, we can rewrite (68) in terms of fi as 

{0])^,T = {O]) [I - 0.009Aefr£^ (/. - /icj + , (71) 

for fi — ficr ^ fJ'cr- Keeping the temperature low while increasing the density, the system stays 
in the same phase. In order to obtain {Oj)d,T in this regime, the integrals in (62) and (64) 
should be evaluated numerically. Figure 7 shows the behavior of the xSB order parameter 
{Oj)d,T/{Oj) as well as Ui — ut at low temperatures but arbitrary density. As was the case for 
the zero temperature, at fixed temperature both of these quantities decrease with d for small 
d and then increase for larger d. On the other hand, at a fixed small density d, {0'^^d,T/ {Oj) 
increases with temperature, whereas Ut — ut decreases with temperature at a fixed density 
of arbitrary value. However, from Figure 7a it appears that (C])d,T/ (C*}) actually decreases 
with temperature for a fixed density at a larger value. Thus, the behavior of 
and Ut — Ut are qualitatively the same (opposite) for smaller (larger) densities. 




Figure 7: (a) {Oj)d,T/ (Cp as a function of T and d in a xSB phase of the holographic NJL model at low 
temperatures but arbitrary density, (b) ut — ut as a function of T and d in the same phase. 



3.3 Finite density and temperature: xSR phase 

At high enough temperatures and densities, parallel branes which support a non-vanishing 
density are thermodynamically favored. Thus, the holographic NJL model at high tempera- 
ture and density is in a xSR phase. Since Oj is charged under chiral symmetry, one should 
be able to show that {0'^^d,T vanishes identically for any contour at m = ua, given that the 
string boundary conditions on parallel flavor branes are preserved. Although it is challenging 
to explicitly show that {Oj)d,T vanishes in general for any contour, in this section we consider 
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a contour that is straight along the x direction and analyze the corresponding worldsheets in 

the xSR phase of the model. The string equation of motion admits a solution for which the 
worldsheet only extends in the u and x directions and is located at fixed points in the rest 
of the directions. Such a worldsheet does not couple to Aq on the flavor branes and is the 
simplest possible worldsheet one can study. Since the flavor branes do not form a U-shaped 
conflguration above the horizon, this worldsheet has to go through the horizon. To better 
analyze the worldsheet, we flrst consider the situation for which the background geometry 
(10) has Minkowski signature. We then consider the worldsheet in the Wick-rotated version 
of the background geometry. This will enable us to compute the thermal expectation value 
of the dual operator in the xSR phase (with flnite density) of the holographic NJL model. 

In order to study the worldsheet bounded by parallel branes in the background geometry 
(10), it is convenient to use Kruskal coordinates since this enables us to see what happens to 

— 1/2 

the worldsheet as it passes through the horizon. First, let's scale u utu and ^ Uj. ' 
in (10) so that the horizon is now located u—1. Next, deflne the Kruskal coordinates by 

^ = +e3(*+0/2^ ^ ^ _e-3(t-.)/2^ (72) 

where r is related to u through 

d^ = ^^^^ 

Note that r ~ | ln(M — 1) near m = 1, so in the new coordinate the horizon is located along 
the V and w axes. The singularity, on the other hand, corresponds to the curve vw — 1. In 
Kruskal coordinates, the background metric takes the form 

ds^ 

— — = —p(v, w)dvdw + q(v, w)dx^ + siv, w^dO^, (74) 
7 

where 7^ = ^TR^, and p{v,w), q{v,w) and s{v,w) are related to the Schwarzchild-like 
coordinate u through 

p{v, w) = ^ "^3^/ , q{v, w) = u^^^, s{v, w) = li^/^ (75) 

Since vw = —e^^, p{v,w), q{v,w) and s{v,w) are all functions of vw only. As we alluded 
to earlier, the simplest worldsheet is extended only in the u and x directions. Suppose that, 
with regards to the Schwarzschild time coordinate, the worldsheet is at t = to- In Kruskal 
coordinates, this translates into w — —e~^*°v. Since the worldsheet does not depend on time, 
one can shift t = to to t = 0, which results in the worldsheet being described by the line 
w = —V in Kruskal coordinates. Thus, the worldsheet, depicted in Figure 8 by the horizontal 
line, passes through the bifurcate horizon and goes to the second asymptotic region in the 
extended geometry where it has nowhere to end. In other words, the worldsheet has inflnite 
area, which results in (Oj) —0. 

The analytic continuation of (Oj) to Euclidean signature should be interpreted as a 
thermal expectation value of the operator. Thus, upon Wick rotation, we should again flnd 
that the worldsheet has inflnite area, which results in a vanishing thermal expectation value 
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singularity 



Figure 8: The Kruskal diagram of the background geometry. The horizontal line depicts the infinite-area 

worldsheet. 

of the dual operator. Let's see this exphcitly. To analyze the behavior of the worldsheet, it 
suffices to just look at the t — u plane, since the x direction simply comes along for the ride. 
Near u = 1, the Euclidean version [dt'^ — > —dt^) of the background metric (10) in the t — u 
plane takes the form 

^^dp' + p'de', (76) 

7 

where p — 2\^u — and 9 — Me/'^- Since the period of is given by (11), the period 

of 9 is 27r and the geometry is smooth ai p — Q {u — 1). In the {p, ^}-coordinates, the 
worldsheet is located at pA in the UV region and is at a fixed value of 6, say ^ = 0. This 
worldsheet passes through the horizon and emerges at the antipodal point 9 = 27r/3, where 
there is no O insertion on which it can end. Equivalently, the worldsheet has infinite area, 
which results in (Cp = 0. This imphes that the xSB order parameter is discontinuous at the 
tri-critical point on the phase diagram of the holographic NJL model at finite temperature 
and chemical potential (Figure 5), where the three distinct phases of xSB with zero density, 
xSB with non-zero density and xSR with non-zero density are equally thermodynamically 
favored. 
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